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Andreev bound states are an expression of quantum 
coherence between particles and holes in hybrid 
structures composed of superconducting and non- 
superconductmg metallic parts. Their spectrum carries 
important information on the nature of the pairing, 
and determines the current in Josephson devices. Here 
I focus on Andreev bound states in systems involving 
superconductors and ferromagnets with strong spin- 
polarization. 1 provide a general framework for non¬ 
local Andreev phenomena in such structures in terms 
of coherence functions, and show how the latter link 
wave-function and Green-function based theories. 


1. Andreev reflection phenomena 

In an isotropic non-magnetic superconductor the normal- 
state single-particle excitation spectrum is modified in 
the superconducting state to — [(e^ — 
acquiring a gap A around the electrochemical potential /r, 
and the density of states is characterized by a coherence 
peak just above the gap, accounting for the missing 
sub-gap states. This spectral signature, predicted by 
Bardeen-Cooper-Schrieffer (BCS) theory [1,2] has been 
first observed by infrared absorption spectroscopy [3,4] 
and by tunneling spectroscopy [5]. 

The spectral features above the gap may show 
information about electron-phonon interaction (or about 
interaction with some low-energy bosonic modes, 
e.g. spin-fluctutations [6]), or may exhibit geometric 
interference patterns. Features due to electron-phonon 
interaction, predicted by Migdal-Eliashberg theory [7, 
8] and studied in detail by Scalapino et al. [9], were 
measured early in tunneling experiments by Giaever et 
al. [10]. They are a consequence of electronic particle-hole 
coherence in a superconductor and build the basis for the 
McMillan-Rowell inversion procedure for determining 
the Eliashberg effective interaction spectrum a^F{u) 
[11]. Geometric interference effects include oscillations 
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in the density of states, such as Tomasch oscillations in N-I-S and N-I-S-N' tunnel structures with 
a superconductor of thickness dg (and with transverse Fermi velocity wf,s)/ giving rise to voltage 
peaks at eVn = [(2Z\)^ + (mrhvp^s/dg)'^]^^^ (n integer) [12-14]; and Rowell-McMillan oscillations 
in N-I-N'-S tunnel structures with a normal metal N' of thickness djqi (transverse Fermi velocity 
uf,n')a giving rise to voltage peaks at eVn = nnhvp^pii/2dpii [15]. Possible offsets due to spatial 
variation of the gap A may occur [16[. 

Inhomogeneous superconducting states exhibit also features at energies inside the gap. Surface 
bound sfates in a normal metal overlayer on a superconductor were predicted first by de Germes 
and Saint-James [17,18] and measured by Rowell [19] and Bellanger et al. [20]. These de Gennes- 
Saint-James bound states have a natural explanation in terms of the so-called Andreev reflection 
process, an extremely fruitful physical picture suggested in 1964 by Andreev [21]. For example, 
geometric resonances above the gap appear due to Andreev reflection at N-S interfaces, which 
describe (for normal impact) scattering of a particle at wavevector kp ^; -F (E"^ — fhvp^x 

into a hole at wavevector kp ^: — {E^ — jhvp^x or vica versa. Below the gap Andreev 

reflections lead to subharmonic gap structure due to multiple Andreev reflections at voltages 
Vn = 2Afen (n integer) in SIS junctions [22] (for a general treatment in diffusive systems see 
[23]), and control electrical and thermal resistance of a superconductor/normal-metal interface 
and the Josephson current in a superconductor/normal-metal/superconductor junction. The 
Andreev mechanism also gives rise to bound states in various other systems with inhomogeneous 
superconducting order parameter, which are named in the general case Andreev bound states. 

Transport trough an N-S contact is strongly influenced by Andreev scattering, and is described 
in the single-channel case by the theory of Blonder, Tinkham, and Klapwijk [24], generalized 
to the multi-channel case by Beenakker [25]. Andreev scattering at N-S interfaces is the cause 
of the superconducting proximity effect [26,27]. Interference effecfs in transport appear also as 
the result of impurity disorder. In contrast to imconventional superconductors, where normal 
impurities are pair breaking, isotropic s-wave superconductors are insensitive to scattering from 
normal impurities for not too high impurity concentration, which is the content of a theorem by 
Abrikosov and Gor'kov [28,29], and by Anderson [30]. In strongly disordered superconductors 
(weak localization regime) the superconducting transition temperature Tc is reduced [31], 
accompagnied by localized tail states (similar to Lifshitz tail states in semiconductors [32]) 
just below the gap edge [33,34]. An interference effect is the so-called reflectionless tunneling 
[25,35-38], which leads to a zero-bias conductance peak in a diffusive N-TS structure. It results 
from multiple scattering of Andreev-reflected coherent particle-hole pairs at impurities, and from 
the resulting backscattering to the interface barrier, making the barrier effectively transparent near 
the electrochemical potential for a pair current even in the tunneling regime. 

Abrikosov and Gor'kov developed in 1960 a theory for pair-breaking by paramagnetic 
impurities, showing that at a critical value for the impurity concentration superconductivity is 
destroyed, and that gapless superconductivity can exist in a narrow region below this critical 
value [39]. Yu [40], Shiba [41], and Rusinov [42] (who happened to work isolated from each other 
in China, Japan, and Russia) independently discovered within the framework of a full t-matrix 
treatment of the problem that local Andreev bound states (now called the Yu-Shiba-Rusinov 
states) are present within the BCS energy gap due to multiple scattering between conduction 
electrons and paramagnetic impurities. Andreev bound states also exist in the cores of vortices in 
type II superconductors. These are called Caroli-de Germes-Matricon bound states [43], and carry 
current in the core region of a vortex [44]. Their dynamics plays a crucial role in the absorption of 
electromagnetic waves [45-47]. 

In an S-N-I or S-N-S jimction, Andreev bound states appear in the normal metal region at 
energies below the gaps of the superconductors. The number and distribution of these bound 
states depend on details such as interface transmission, mean free path, and length of the normal 
metal dpi. In general, there is a characteristic energy, the Thouless energy [48] (related to the dwell 
time between Andreev reflections), given by hvp pi/dpi for the clean limit, and by hDpi/d'^ for 
the diffusive limit, with Fermi velocity wf,n and diffusion consfant 79^ of the normal metal. 
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In the diffusve limit, the Andreev states build a quasi-continuum below the superconducting 
gap, whereas in the case of ballistic junctions bands of Andreev bound states arise. For the case 
that no superconducting phase gradient is present in the system, however, a low-energy gap 
always arises in the spectrum of Andreev states in the normal metal. This so-called minigap 
scales for sufficiently thick normal metal layers approximately with its Thouless energy and 
with the transmission probability (possibly further reduced by inelastic scattering processes). It 
was found first by McMillan [49] and can be probed by scanning tunneling microscopy [50]. In 
chaotic Andreev billiards [51], where disorder is restricted to boundaries, a second time scale, the 
Ehrenfest time, competes with the dwell time to set the minigap [52]. 

The importance of Andreev bound states in S-N-S Josephson junctions for current transport 
was first discussed by Kulik in 1969 [53]. Andreev bound states form in a sufficiently long normal 
region, which are doubly degenerate (carrying current in opposite direction) for zero phase 
difference between the superconducting banks. For a finite phase difference, this degeneracy is 
lifted. The gap or minigap in a Josephson structure is reduced and eventually closes when a 
supercurrent flows across the junction. This is a result of a "dispersion" of fhe energy ii'b.s. of the 
Andreev bound states as function of phase difference Ax between the superconducting banks 
[53-55]. The contribution of the bound states to the supercurrent is given by {^el'fi)dE\, ^ jdAx, 
with e < 0 the electron charge. Apart from the current carried by the Andreev bound states, there 
is also a contribution from continuum states above the gap [53]. For a single-channel weak link 
between two superconductors with normal-state transmission probability there is one pair of 
Andreev boimd states with dispersion E± = ±Zl[l — sin^(Zix/2)]^/^. 

The large size of a Cooper pair in conventional superconductors leads to a pronoimced 
non-locality of Andreev reflection processes. This allows for interference effects due to crossed 
Andreev reflection, in which the particle and hole involved in the process enter different normal- 
state (typically spin-polarized) terminals, which are both simultaneously accessable to one 
Cooper pair [56,57]. This effect has been first experimentally observed by Beckmann et al. [58]. 

Finally, an important role is played by Andreev zero modes as topological surface states. 
Examples are zero-bias states at the surface of a d-wave superconductor [59,60], and Majorana 
zero modes in topological superconductors [61,62] and superfluids [63]. 


2. Andreev bound states at magnetically active interfaces 


(a) Spin-dependent interface scattering phase shifts 


The importance of spin-dependent interface scattering phase shifts for superconducting 
phenomena has been pioneered in the work of Tokuyasu, Sauls, and Rainer in 1988 [64]. Consider 
an interface between a normal metal (N) at a; < 0 and a ferromagnetic insulator (FI) or a half- 
metallic ferromagnet (HM) at a; > 0. For simplicity, let us model the FI (or HM) by a single 
electronic band with energy gap for spin-down particles and an energy gap = V^ — 2 J 
for spin-up particles, where J > 0 denotes an effective exchange field. The exchange field can 
be related to an effective magnetic field via fjSeS = J (for free electrons the magnetic moment 
is fi = fie <0). Let us assume an incoming Bloch electron with energy 0 < E <Vf and spin 
(TGlt,),}, reflected back from the interface with amplitude ra- It is described by a wave 
function Ecrix, ry) = (e*^“ -I- To- e“'^®) at a; < 0 and Eaix, ry) = at a; > 0. For 

the normal metal hk{E) = [2mE — (fifey)^]^/^. For the FI hKa{E) = [2m(lo- — S) -I- (fifcy)^]^/^. 
The reflection scattering matrix is 


S = 


A'&'f 


0 



^ k + ’ 



k — 
k -I- 


( 2 . 1 ) 


In the range < E <Vj^ the spin-up electron can be transmitted for sufficiently small fcy with 
amplitude = 2^kk^/{k -\- k^), where fifc^(i5) = [2m(i? — V^) — (fifcy)^]^/^. In this case, the 
reflection amplitude is also real, and equal to Va = {k — k^)/{k + k^). The reflection phase is —tt 
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Figure 1. (a) Spin dependent scattering phase shifts for Bloch waves with energy Ep reflected from an N-FI or N-HM 
interface (at x = 0). Here Vi = 3Ep is fixed and varied from —Ep to 3 -Bf, fc|| = 0. For Vf > Ep this describes 
an FI, for Vf < Ep a HM. For Vff < 0 the Fermi surface of the spin-up band in the HM becomes larger than in N. 
Shown are for x < 0 the (normalized) reflected wave, quantified by —lm(rc e“''“^)/|ro-|, and for a; > 0 the transmitted 
wave, quantified by —(for real Ko-) or —(for real k-^, i.e. Vf < Ep). The spin-up wave is 
shown in orange, the spin-down wave in blue, (b) Spin mixing angle •& as function of parallel momentum fc||/A:F and 
u = {Vf + Vf)/2 for J = O.SEf (FI for v > 1.3, HM for v > 0.7); (o) the same for J = O.SEf (FI for v > 1.8, HM for 

IT >0.2). 


for k <k^ and zero for fc > In figure 1 (a) = ZE-p is fixed and V-f varied from —Ep to ZEp, 

fc|| = 0. A phase shift for reflected waves between the two spin projections appears. It results 
from the well-known effect that reflection from an insulating region results in a phase-delay of 
the reflected wave with respect to the case of an infinite interface potential. This phase delay 
appears due to the quantum mechanical penetration of the wave function into the classically 
forbidden region. The range V-f — Ep > 0 in figure 1 (a) corresponds to an N-FI interface, where 
both spin-projections are evanescent in FI. Here, the reflected spin-up wave trails that of the spin- 
down wave, and the effect increases when — Ep approaches zero. The phase 'd = '&^ — of 
the parameter called spin-mixing angle [64], or spin-dependent scattering phase 

shift [65]. It is an important parameter for superconducting spintronics. For the N-FI model 
interface it is given by 


d d^ — d\ 

tan — = tan —-^ 

2 2 


Hh - «t) 

_|_ Kf-Kff ’ 


( 2 . 2 ) 


which is positive due to The range V-f — Ep < 0 corresponds to a N-HM interface, with 

the spin-up band itinerant in HM. Here, as long as the spin-up Fermi surface in the HM is smaller 
than that in N (for — 1 < (H^ — Ep)/Ep < 0), the reflection phase in is zero, and 


-d —d\ K 

tan - = tan-k = ^ (2.3) 

2 2 fc ’ ^ ^ 

which can aquire large values. Finally, in the range {V^ — Ep)/Ep < —1 the Fermi surface in the 

HM is larger than in N, which leads to a reflection phase of vr for spin-up particles, and 


d 

tan — = tan 



A 


(2.4) 


which now is negative. In ballistic structures, the spin-mixing angle depends on the momentum 
Iik|| parallel to the interface, as illustrated in figure 1 (b) and (c) for varying Fermi surface 
geometry in the ferromagnet, here parameterized by varying {V^ -F V\f)/2 keeping J fixed. If both 
spin-bands are itinerant in the ferromagnet (F), then the spin-mixing angle is zero (if fe > or 

k < kf-, kjf) or —TV (if k^ > k> k\^, red areas in figure 1 b and c), unless an interface potential exists, 
rendering the reflection amplitudes complex valued. In general, the spin-mixing angle should be 
considered as material parameter, which in addition depends on the impact angle of the incoming 
electron or on transport channel indices. 

Note that the parameter r^r* has become well-known in the spintronics community, as it 
governs the spin mixing conductance [66] in spintronics devices. 


rsta.royalsocietypublishing.org Phii. Trans. R. Soc. A 0000000 






5 


It is also instructive to study an incoming Bloch-electron polarized in a direction different 
from the magnetization direction in the ferromagnet. Let us consider the case of a FI. For a Bloch 
electron polarized in a direction n(Q;, (ji), parameterized by polar and azimuthal angles, a and (j), 


\ * 


-ikii 


ikx 


COS 



tz + sin 



J-zj 


gikyriigifea: 


(2.5) 


the reflected wave will have the form 

r a i, . a 

cos—e 2 t^+sin—e ^ 


iz] e‘^ 


— ikii 

e II 


I —\kx _X ii? ikii 

'e e " 


— ikx 


( 2 . 6 ) 


with i? = + i?^)/2. Similarly, 4,^ ^ scatters into This means that scattering leads, 

apart from an unimportant spin-independent phase factor to a precession of the spin around 
the magnetization axis [64]. The direction of precession depends on the Fermi surface geometries, 
and is determined by the sign of the spin-mixing angle 


(b) Andreev reflection in an S-N-FI structure 

An important consequence of spin-mixing phases is the appearence of Andreev bound states at 
magnetically active interfaces, predicted theoretically [67-72,77], and verified experimentally [78]. 

Consider a superconductor near an interface with a ferromagnetic insulator. Let us assume 
that the superconducting order parameter is suppressed to zero in a layer of thickness d next to 
the FI interface, such that the structure can be described as a S-N-FI junction with identical normal 
state parameters in S and N. For simplicity I consider here a spatially constant order parameter in 
S (extending to the half space a: < 0). The FI (a; > d) will be parameterized by reflection phases 
andwith spin-mixing angle Solving the corresponding Bogoliubov-de Germes 

equations in the superconductor (ui are spin Pauli matrices, (jq a 2 x 2 unit spin matrix) 




2m 

-A 


icr 2 


Aicr 2 



(2.7) 


with spinors u and v, the (still unnormalized) eigenvectors for given energy e and ky = 0 are 




( M 




7 _7 \ 

0 

6 ^ , 

1 

■^ikj^x 

e ^ , 

7 

±i/c_ X 

e , 

0 

0 


-7 


1 


0 

V 7) 




(0) 


( 1 / 


( 2 . 8 ) 


where to first order in \A\/E-p the wavevectors are k±{e) = kp ± i-y/|Z\p — j(hvp) and 

A ,, , A* 


7(e) = - 


e + VlAp -e 2 ’ 


7(e) = + 


e + i-v/|A |2 - £2 


(2.9) 


For I e I > IAI the expression i Ap — ^2 jg replaced by e \J\ — (| A|/e)2 (which corresponds to e —> 
£ -F iO'*' with infinitesimally small positive 0^). In the N layer the solutions are obtained by setting 
A = 0. In the FI only evanescent solutions of the form and are allowed. The reflection 

coefficients connect incoming (e**+^, solutions with outgoing (e“'^+“, e‘^““) solutions. 

For scattering from electron-like to electron-like Bogoliubov quasiparticles and for electron-like 
to hole-like Bogoliubov quasiparticles in leading order in \A\/Ep they are 


ret-i-et ^ 




4idT 


1 + 776 e 


id 


g-2idfcFg2i'iR%g-idr 

rh'[^h-\ = - _ 4 idpC” 

1 -F 776 '‘”6 6 


(1 + 77) 




— /-I — i'i9\ 

— 7(1 — e e ) 

- , - :,Q ’ 

\ ^ nv-p git7 

— 7(1 — e e ) 

- , - 4id^:^ _:,Q ’ 

\ + 776 '"^F e 


( 2 . 10 ) 

( 2 . 11 ) 


and similar relations hold for and simultaneously d —>■ —d, 7 —>■ — 7 , 7 —>• — 7 . These relations 

have a simple interpretation. The coherence functions 7 and 7 represent probability amplitudes for 
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Figure 2. Imaginary part of Andreev reflection amplitudes for spin-up Bogoliubov quasiparticle to spin-down Bogoliubov 
quasihole, lm(re-t'^r,^), at normal impact, for an S-N-FI structure with a normal region of thickness d, as function of 
energy e, and of d (in units of go = with the projection of the Fermi velocity on the surface normal). The 

reflection amplitude at the N-FI interface is for spin-up e‘’^t and for spin-down The spin-mixing angle is defined 
as = 1 ?^ — 1 ?^. It has the values (a) = 0, (b) = p/2, (c) = tt. In (d)-(f) the thickness of the normal layer is fixed 

to (d) d = 0, (e) d = go. (f) d = 2go, and vary the spin-mixing angle d. The negative reflection amplitude for spin-down 
quasiparticle to spin-up quasihole, —lm(rej,_>(j^), is obtained by inverting the energy axis, e -n- —e. 


id—^ 

hole-to-particle conversion (— 7 ) or particle-to-hole conversion (— 7 ), whereas the factors e '“"f 
represent the electron-hole dephasing when crossing the N layer. Thus, the factors 776 ^'“^ '‘"f 
represent a Rowell-McMillan process of four times crossing N with two reflections from FI (once 
as particle and once as hole, contributing and and two Andreev conversions. When 

this factor equals —1, which happens for energies below the gap, a bound state appears in N 
due to constructive interference between particles and holes. Note that for |£| < |Z\| the coherence 
functions have unit modulus: I 7 I = I 7 I = 1, such that with A = |Zi|e‘^ one can write 7 = 
and 7 = —withsin!!/= e/Izi|, cosll/> 0. For |e| < |zi| only outgoing wavevectors—fe+ 
and k- lead to normalizable solutions in the superconductor, which are restricted to the bound 
state energies, given by the solution of e = | Aj cos(|^ ± f )sign[sin(|^ =•= f )]• 

In figure 2 results for the quantity Im(rg^_^^^) for normal impact are shown. For energies 
above the gap the typical Rowell-McMillan oscillations are visible. Below the gap sharp bound 
states exist, the energy of which depends on bofh and d. In (a)-(c) the influence of varying d is 
illustrated. The special case i? = 0, shown in (a), corresponds to the classical Rowell-McMillan 
S-N-I structure. For = tt, shown in (c), a midgap boimd state exists for all d. For 0 < i9 < tt 
particle-hole symmetry is broken, as seen in figure 2 (b) and (d-f). A corresponding bound sfate 
at negative energy exists for The two corresponding bound states in the density of states 

have opposite spin polarization. With increasing d more and more bound states enter the sub¬ 
gap region, emerging from the continuum Rowell-McMillan resonances. For d=0 a bound state 
exists for any nonzero i9, as shown in (d), and as discussed in Ref. [67]. The influence of is 
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shown for various d in (d)-(f). These results can be interpreted as one spin-polarised chiral branch 
crossing the gap region with increasing For d = 0 this happenes when varying d from zero to 
27r. For d / 0 one needs a variation exceeding 27r (up to multiple times) until the branch crosses 
the entire gap. The figure shows results for normal impact, fey = 0. In general, an integration 
over fcy will lead to Andreev bands instead of sharp bound states, similar as in the case of de 
Gennes-Saint-James bound states in S-N-1 structures. 

Finally, note that with the reflection matrix (2.1) the resulting coherence function developes a 
spin-triplet component from a singlet component 7 '" = 7 oi(T 2 : 


( ° 
I 



( 2 . 12 ) 



which implies that a singlet pair is scattered into a superposition of a singlet and a triplet pair: 


(t; - it) ^ (t; e*’’- it e-*’’) = cos(d)(ti - it) + i sm(d)(ti + it)- (2.13) 


(c) Andreev bound states in an S-FI-N structure 

As next example 1 summarize some results from Refs. [79,80] and section IV of Ref. [82] for an 
S-Fl-N junction, consisting of a bulk superconductor coupled via a thin ferromagnetic insulator 
(such as EuO) of thickness di to a normal layer of thickness . 1 assume here the ballistic case, 
and refer for the diffusive case to Refs. [79-81]. The interface is characterized by potentials 
and Vi = Vf + 2J, such that the energy dispersions in the superconductor is I 2 ra, in the 
normal metal Vn + j‘Ira, and in the barrier Va -F j^m, a £ {t, 4 -}- The parameter Vn is 
used to vary the Fermi surface mismatch. The Fermi wave vectors and Fermi velocities in S and 
N are kp s, vp 3 and kp n, vp n, respectively. The Fermi energy is E-p = s/ 2 ™- For shorter 

notation we introduce a directional vector for electrons moving in positive i-direction, vp n = 
|vp N|h (i.e. fix > 0), and the corresponding i-component of the Fermi velocity, uf_n,x = Vx> 0 , 
in the normal metal (situated at 0 < a; < djsf). It is convenient to define coherence functions 7 and 
7 as 2x2 spin matrices. The coherence functions in the superconductor are given by 7 = 7oio'2 
and 7 = 7oicr2, where 70 and 70 are given by the expressions in (2.9). The solutions in the normal 
metal are (for simplicity of notation 1 also suppress the arguments ky and e in 7 and 7) 

'rinx,x) = 'y{hx, , 'y{-hx,x) = 'y{-nx, (2.14) 

7(-ha;,a:) =7(-hz,,0)e^‘^'^/^"", 7(ha:, x) = 7(nx, (2.15) 


At a; = dN one obtains 7(h£c, dN) = 7(— ha;, dN) = 7 b arid ^{fix, djvj) = ^(—fix, djq) = 7b, with 



(2.16) 


The scattering parameters are the modulus of the transmission amplitudes, and the modulus 
of the reflection amplitudes = (1 — and = (1 — (equal on both sides of the FI), 

as well as the phase factors of the scattering parameters (all these parameters depend on ky). The 
relevant energy scale in the normal metal for given direction h(ky) is 


d(ky) = t 7 (ky) ^^(ky) ha;(ky) £Th, ^Th = hWF.N/SdN, (2.17) 

with the Thouless energy exh- Matching the wavefunctions at x = 0 to the thin FI layer and the 
superconductor, leads to 7 - = —7-(- and 7 + = — 7 -, as well as [79,80] 

5 


= - 


Va + — {Va + i0+)2 


(2.18) 
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Figure 3. Energy-resolved DOS in the normal metal for different values of the interface exchange field J. The energy 
scale is So = (t-t-*J.^Th)k||=o. the Thouless energy exh = ftrrp.N/SdN- The interlayer thickness is d\ = 2/kp^s 
and the interface potentials are = 1.2Ep, Vi = V^ + 2J. The width of the normal layer is ciN = j A. The inset 
in the lower left corner of each panel illustrates the Fermi-surface mismatch: in (a) = O.Sfcp s, in (b) fcp,N = ^f.s. 

and in (c) fcp^N = lOfcp^s- Adapted from [80], Copyright (2010) by the American Physical Society. 


where cr G {+i aird the function Va(s) is defined as 


txo-(e) = hxSTh 


riajETh 


+ crd+ + < 1 ' 


+ sin 


VWa;£Th 


+ crd- 



(2.19) 


with 'd± = where i?n arid i9g are the spin-mixing angles for reflection at the Fl-N 

interface and the S-FI interface, respectively, and the variable a is to be understood as a factor ±1 
for (J = ±. Note that (2.18) has the same form as (2.9) with the role of A and e taken over by 5 and 
Vcr, respectively. Note also that | 7 o-| = 1 for Ucr < S, even in fhe funneling limit. This is an example 
of reflectionless tunneling at low energies and resuls from multiple reflections within the normal 
layer. Quasiparticles in the normal layer stay fully coherent in this energy range. 

The density of states at the outer surface of the N layer is obtained as 


NBje) 

Nf,N 




l + 7a 
i-7i 


= Re^ 


\Va{e)\ 




( 2 . 20 ) 


where {...) denotes Fermi surface averaging, and Np n is the density of states at the Fermi level of 
the bulk normal metal. Results for this density of states are shown in figure 3. The various panels 
(a)-(c) show examples for various Fermi surface mismatches. In (c) there are non-transmissive 
channels present in the normal layer (|kj| | > kp s), leading to a large constant background density 
of sfafes. In each panel, the curve for J — 0 corresponds to the case of an non-spinpolarized SIN 
junction. There is a critical value Jcrit (independent from the Fermi surface mismatch and equal 
to « O-lSSp in the figure) above which fhe system is in a state where no singlet correlations are 
present in the normal metal at the chemical potential (e = 0), and pure odd-frequency spin-triplet 
correlations remain. In this range the density of stafes is enhanced above its bulk normal state 
value [80]. On either side of this critical value the density of states decreases as function of J, 
however stays always above Np n for J > Jcrit- In the diffusive limif a similar scenario arises, 
with a peak centered at zero energy in the density of states [79,80]. A zero-energy peak in the 
density of states has been suggested as a signature of odd-frequency spin-triplet pairing also 
in hybrid structures with an itinerant ferromagnet or a half-metallic ferromagnet coupled to a 
superconductor [83-85]. 

It is interesting to study the tunneling limit, ^ 1, <C 1, for small excifation energies e ^ 
min(extu 2i) and small spin-mixing angles Jn/ '^S- Then Ua ~2£ + crfixeThi^N/ i-e. i^a depends 
only on fhe spin mixing angle af the Fl-N interface, which acts in this case as an (anisotropic) 
effective exchange field b = ha;exhi?N/2 on the quasiparticles. For diffusive sfructures, a similar 
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connection between an effective exchange field and the spin-mixing angle has been made [ 86 ]. 
The parameter 5/2 = on the other hand acts as effective (anisotropic) gap function. 

For each direction n, the gap closes at a critical value of effective exchange field, 6 = 5/2, which 
happens for • 

3. Andreev bound states in Josephson junctions with strongly 
spin-polarized ferromagnets 

(a) Triplet rotation 

Interfaces with strongly spin-polarized ferromagnets polarize the superconductor in proximity 
with it, as shown in the previous section. However, inorder for superconducting correlations to 
penetrate the ferromagnet, it is necessary to turn the triplet correlations of the form fi + it info 
equal spin pair correlations of fhe form ft and ii. The reason is that correlations involving spin- 
up and spin-down electrons involve a phase factor fcpt “ which in strongly spin-polarized 
ferromagnets oscillates on a short length scale. This leads to destructive interference and allows 
fo neglect such pair correlations on the superconducting coherence length scale [87]. 

The way to achieve this is to allow for a non-trivial magnetization profile at the interface 
between the ferromagnet and the superconductor. This can include for example strong spin- 
orbit coupling, or a misaligned (with respect to the bulk magnetization) magnetic moment in 
the interface region. For strongly spin-polarized ferromagnets this has been suggested in Ref. 
[68,88,89]. For weakly spin-polarized ferromagnets a theory was developed in 2001 involving a 
spiral inhomogeneity on the scale of the superconducting coherence length [90,91]. A multilayer 
arrangement was subsequently also suggested [92,93]. For various reviews of this field see Refs. 
[87,94-104]. 

The idea is to rotate the triplet component, once created by spin-mixing phases in the S- 
F interfaces, into equal-spin triplet amplitudes with respect to the bulk magnetization of the 
ferromagnet [101]. This is achieved by writing a triplet component with respect to a new axis 

(t/+ = ~ sin(a) l^e '^(tt)^ ~ e +cos(Q)(t4.-l-4.t)^i (3-1) 

where a and (j> are polar and azimuthal angles of the new quantization axis. Then, if a thin FI layer 
oriented along the {a, <p) direction is inserted between the superconductor and the strongly spin- 
polarized ferromagnet with magnetization in z direction, equal spin-correlations can penetrate 
with amplitudes — sin(Q)e“*'^ and sin(a)e“^, respectively. These correlations are long-range and 
not affected by dephasing on the short length scale associated with — kpj^. 

(b) Pair amplitudes at an S-FI-F interface 

It is instructive to consider the scattering matrix of an S-Fl-F interface between a superconductor 
and an itinerant ferromagnet, with a thin FI interlayer of width d, in the tunneling limit. In this 
case one can achieve an intuitive understanding of fhe various spin-mixing phases involved in 
the reflection and transmission processes. Denoting wavevector components perpendicular to 
the interface as k in the superconductor, and in the ferromagnet (1 assume ky such that both 
spin directions are itinerant), and imaginary wavevectors and in the FI, the FI magnetic 
moment aligned in direction (sin a cosy;, sin a simp, cos a), and a F magnetization aligned with 
the 2 -direction in spin space (a = 0), matching of wavefunctions leads to a scattering matrix 

s. { . j f 1, \ { , j (3.2) 

\ 0 J V / V 0 J 

where ^’s.F are phase matrices which include the spin-mixing phase factors. Dip, Da, Dp 
are spin-rotation matrices, i's,F carry information about S-Fl and Fl-F wavevector mismatch. 
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and T contains the tunneling amplitudes including wavevector mismatch between S and F. 
In particular, if one denotes diagonal matrices with diagonal elements a, b by diag(a, b), then 
K — diag{K^/k, K^/k), Q = diag{q^/k,q^/k), the spin-rotation matrices Dip, Da between the 
quantization axis in the FI and the 2 ; axis, and the phase matrices <^s,F are 


D 


ip — 




cos ^ — sin ^ 

sin ^ cos § 


, ^S.F = 




0 


0 ^ 

) ’ 


(3.3) 


and the spin-rotation matrix Dp at the FI-F interface results from Q~^DaK dI,Q~^ = DpZD^^ 
with Z = diag(^^, The angle fi vanishes for a = 0, and varies from K^/q^ at a = 0 to 
Kp/q^ at a = TT, correspondingly varies from K,p/qp to n^/qp. Also, = (1 — iiT)/(l-I-iiT), 
^F = (1 — + '^Z), up = ^2I{1 + K'^), Up = .y/ 2 /(l -|- ^ 2 ), and the tunneling amplitude is 

T =V A with V = diag(e“''^‘^, and the real-valued mismatch matrix A — KDaQ~ 2 Dp = 

DaQ^DpZ, oi which the off-diagonal elements appear for a 7 ^ 0, tt only. 

One can see from equation (3.2) that the spin-mixing phases, which appear in the reflection 
amplitudes, also enter the transmission amplitudes; in the tunneling limit they contribute from 
each side of the interface one half [ 68 ]. Furthermore, one should notice that the interface is 
described by two spin-rotation matrices: one given by the misalignment of the FI magnetic 
moment with the 2 axis in spin space, and one which is combined from the magnetization in F and 
the magnetic moment in FI. The latter appears because the wave function at the FI-F interface is 
delocalized over the FI-F interface region on the scale of the Fermi wavelength and experiences an 
averaged effective exchange field, which lies in the plane spanned by the 2 axis and the direction 
of the FI magnetic moment (same Dip in equation (3.2)). 

Pair correlation functions / are related to coherence fimctions by / = —27ri7(l — 77 )”^, with 
/, 7 and j 2x2 matrices in spin space [105]. When both are small (near Tc or induced from a 
reservoir by tunneling through a barrier), / and 7 are proportional. Assuming an incoming singlet 
coherence function 70 in S, the coherence functions reflected back into S and the ones transmitted 
to F can be calculated to linear order in the pair tunneling amplitude according to 


7‘ui^Sn 


0 

-70 


70 

0 


a* (F) Q 

’ 7 out ^21 


0 

-70 


70 

0 


For the reflected amplitude in S one obtains (i9s = 117 “ ) 


(S) / _■ 

7qjj( I —i sin i9g sin ae cos ilg + i cos a sin ilg 

7o I — cos 'dp -F i cos a sin dp i sin dp sin 


(3.4) 


(3.5) 


which is just equation (2.12) rotated in spin-space by the spherical angles a and (p. For the equal- 
spin coherence functions (or pair amplitudes) in F follows up to leading order in the misalignent 
angles a, fi (denoting ilp = ^7 “ 


7 outtt 

70 


l"F7 


(F) 

7out71 


70 


Up^ 


I— sin(a) — sm(/i) 
sin(Q:) — sin(/ 3 ) 


+ sin 


+ sin 


’^S+’^F^sin(/3) 




(3.6) 

(3.7) 


with C = 4 e~^'‘^~'~'^'‘-^'^(up^uppup^upp)K^Kp/[k(qpqp) 2 ]. The transmitted fi and coherence 
functions 7 out 77 ~ and 7out7t ~ spatially oscillate with a 

wavevector 6**^1711 in F, and are suppressed (except in ballistic one-dimensional 

channels) due to dephasing after a short distance 1/| 97 — <? 7 1 away from the interface. Importantly, 
from (3.6) and (3.7) if is visible thaf the equal-spin amplitudes acquire phases ±<p from the 
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azimuthal angle in spin space, which play an important role in Josephson structures with half- 
metallic ferromagnets [87,89] and with strongly spin-polarized ferromagnets when two interfaces 
wifh different azimuthal angles ipi and are involved [106,107]. 

The misalignment of FI with F also induces a spin-flip term during reflection on the 
ferromagnetic side of the interface, which creates for an in F incoming amplitude 7 ^^ a reflected 
amplitude 7 ^^ = sin^ ^ sin^ P, and for an incoming amplitude 7 ^^ a reflected amplitude 

7 |.|. = sin^ ^ sin^ p. In this case, twice the azimutal angle ±ip enters. 

(c) Andreev bound states in S-FI-HM-FI’-S junctions 

For an S-FI-HM interface with a half-metallic ferromagnet (HM) in which one spin-band (e.g 
spin-down) is insulating and the other itinerant, equation (3.6) is modified to [89] 

(F) 

Touttt -IM -{Kf+Ki)d 2 ■ f■ r \ /O QX 

~-4ie ^ — J sin(a). (3.8) 

The same equation also holds at an S-FI-F interface when the conserved wavevector ky is such that 
only one spin-projection on the magnetization axis is itinerant in F. For strongly spin-polarized F 
this is an appeciable contribution to the transmitted pair correlations. 

In order to describe Josephson structures, it is necessary to handle the spatial variation (and 
possibly phase dynamics) of bofh coherence amplitudes and superconducting order parameter 
(which are coupled to each other). For this, a powerful generalization of the 2x2 spin-matrix 
coherence functions has been introduced [82,105], based on previous work for spin-scalar 
functions in equilibrium [108,109] and nonequilibrium [45]. The coherence functions 7 (n, R, e, t) 
fulfil transport equations 

ihvp- V 7 + 2£7 = 7 o o 7 -|- ^17 07 — 70 — A, (3.9) 

\hvp- V 7 — 2ep = p o Ao p -\- 07 — 70 — A (3.10) 

with n a unit vector in direction of vp, X and A include parficle-hole diagonal and off-diagonal 
self-energies and mean fields (e.g for impurity scattering, A includes the superconducting order 
parameter), and external potentials. The time-dependent case is included by the convolution over 
the internal energy-time variables in Wigner coordinate representation, 

(Ao (3.11) 

In the time-independent case it reduces to a simple spin-matrix product. Furthermore, the 
particle-hole conjugation operation is defined by ^(n, R, e, f) = ^*(—n, R, —e, f). Particle-hole 
diagonal [g = —i7r(2V — erg) and g — i 7 r( 2 V — erg)] and off-diagonal [/ = —27riX and / = 27ri.X] 
pair correlation functions (quasiclassical propagators) are obtained in terms of these coherence 
functions by solving the following algebraic (or in the time-dependent case, differential) equations 

V = crg-|- 7 o 7 oV, V = crg-|- 7 o 7 oV, X = 7 -|- 7 o 7 oX, .X = 7 -|- 7 o 7 o.X. (3.12) 

In figure 4 an example for a fully self-consistent calcuation of the spectrum of subgap Andreev 
sfates in a S-FI-HM-FT-S junction is shown, obtained by solving equations (3.9), (3.10) as well 
as the self-consistency equation for the superconducting order parameter in S. For details of the 
calcuation and parameters see [68,96]. The ferromagnetic insulating barriers FI and FT are taken 
identical in this calculation, and the spectra are shown at the half-metallic side of the S-FI-HM 
interface. The mosf prominent feature in these spectra is an Andreev quasiparticle band centered 
at zero energy [68,82,110]. Further bands at higher excitation energies are separated by gaps. The 
zero-energy band is a characteristic feature of the nature of superconducting correlations in the 
half mefal: they are spin-triplet with a phase shift of ±7r/2 with respect to the singlet correlations 
they are created from. The dispersion of the Andreev states with applied phase difference between 
the two S banks determine the direction of currenf carried by these states. This direction is 
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Figure 4. Local density of Andreev states (LDOS) at the HM side of the S-FI-HM interface in a S-FI-HM-FI’-S Josephson 
structure, as function of phase difference Ax over the junction, for quasiparticles with normal impact, at T = 0.05Tc. 
All states are fully spin-polarized, (a) and (c) Dispersion of the maxima of the LDOS as function of phase difference. 
Regions with low LDOS are white, regions of high LDOS (bands of Andreev bound states) are shaded. The signs indicate 
the direction of the current carried by the Andreev states, (b) and (d) show spectra for a fixed phase difference, both for 
positive (full lines) and negative (dashed lines) propagation direction, (a)-(b) is for a large misalignment between FI and 
FIM, and (c)-(d) for a small misalignment, (b) and (d) from [68]. Copyright (2003) by the American Physical Society. 


indicated in the figure by + and — signs. In panels (b) and (d) for a selected phase difference the 
spectra for positive and negative propagation direction are shown. These spectra, multiplied with 
the equilibrium distribution function (Fermi function) determine the contribution of the Andreev 
bound states to the Josephson current in the system. These spectra are shown for normal impact 
direction (ky = 0). An integration over ky gives the local density of states. 

For the case that one can neglect the variation of the order parameter A in S, one can derive 
quite a number of analytical expressions [82,87,111]. Examples for integrated spectra are shown 
in figure 5, taken from Ref. [82]. In (a) the well-known spectrum of de Gennes-Saint James bound 
states is seen for an S-N-S junction [112], showing a dispersion with phase bias Ax between the 
two superconductors. At Ax = tt a zero energy bound state is present, which is a topological 
feature of the particular Andreev differential equations describing this system, for real-valued 
order parameters that change sign when going from the left S reservoir to the right S reservoir. 
The origin is the same as for the midgap state in polyacetylene [113], which is governed by similar 
differencial equations. Such midgap states have been studied in more general context by Jackiw 



Figure 5. Local density of states in the center of a current biased high-transmissive symmetric Josephson junction for 
(a) an S-N-S junction and (b)-(d) an S-FI-FIM-FI-S junction. In (a) and (b) the phase difference Ax over the junction is 
varied from 0 to n. In (b) and (c) the length of the junction is varied for a zero-junction and a tt -junction. The temperature 
is T = O.lTc, the coherence length of the half metal Jo = fi\vp\/2nTc. The FI misalignment angle is a = n/2. The 
transmission parameter t and spin-mixing angle ils depend on the impact angle 'I'n measured from the surface normal; 
this is modeled here by f(<Zrn) = io cos <kn/(l — fg ^ and = jIq cosiTn. Adapted from [82]. Copyright 

(2009) by the American Physical Society. 
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and Rebbi [114] and ultimately have their deep mathematical foimdation in the Atiyah-Patodi- 
Singer index theorem [115]. For the S-Fl-HM-Fl-S junction, shown in (b)-(d), the prominent feature 
for all values of Ax is the band of Andreev states centered around zero energy. The width W 
of this low-energy Andreev band depends on the parameter P = sin (i9s/2) sin(a) and can be 
calculated for fhe limit of short junctions (L —> 0) for f = 1 as [82] 

lR(Ax = 0) = 2|Z\|\/l-P2, W{Ax = tt) = |A|(\/2 - P2 - P). (3.13) 

In the limit P —>• 0 this gives W{Ax = 0) = 2|Zi| and W{Ax = tt) = v^|zi|. Note that compared to 
the S-N-S junction, the low-energy features disperse in opposite direction when increasing Ax for 
fhe S-Fl-HM-Fl-S junction. This means thaf the current flows in opposite direction, and typically 
a TT-junction is realised for identical interfaces. If the azimuthal interface misalignment angles 
(p differ by tt in FI and FT, then this phase would add to Ax according to equation (3.8) and a 
zero-junction would be realized. In the general case, a ((>-junction appears, both for ballistic and 
diffusive structures [87,89,116]. 


(d) Spin torque in S-FI-N-FI’-S’ structures 

Andreev states play also an important role in the non-equilibrium spin torque and in the spin- 
transfer torque in S-F structures [117-120]. Zhao and Sauls found that in the ballistic limit the 
equilibrium torque is related to the spectrum of spin-polarized Andreev bound stales, while 
the ac component, for small bias voltages, is determined by the nearly adiabatic dynamics of 
the Andreev bound states [121,122]. The equilibrium spin-transfer torque req in an S-Fl-N-FT-S' 
structure is related to the Josephson current P, the phase difference between S and S', Ax, and 
the angle Aa between FI and FT, by [123] 


fi^=2e^ 
dAa dAx' 


(3.14) 


Similarly, as the dispersion of the Andreev bound states with superconducting phase difference 
Ax yields fhe confribufion of fhe bound state to the Josephson current, the dispersion of 
fhe Andreev sfates with Aa yields the contribution of this state to the spin current for spin 
polarisation in direction of the spin torque. The dc spin current shows subharmonic gap structure 
due to multiple Andreev reflections (MAR), similar as for the charge current in voltage biased 
Josephson junctions [124,125]. For high transmission junctions the main contribution to the dc 
spin current comes from consecutive spin rotations according to equation (2.6) when electrons 
and holes undergo MAR [122]. 

Turning to ac effects, for a voltage eV <C A the time evolution of spin-transfer torque is 
governed by the nearly adiabatic dynamics of the Andreev bound states. However, the dynamics 
of the bound state spectrum leads to non-equilibrium population of fhe Andreev bound sfates, for 
which reason the spin-tranfer torque does not assume its instantaneous equilibrium value [122]. 
For the occupation to change, the bound state energy must evolve in time to the continuum gap 
edges, where it can rapidly equilibrate with the quasiparticles, similar as in the adiabatic limit of 
ac Josephson junctions [126]. 

The effect of rough interfaces and of spin-flip scattering on spin-transfer torque in the presence 
of Andreev reflections has been discussed by Wang, Tang, and Xia [127]. For diffusive strucfures 
see [128]. Magnetization dynamics has been also addressed recently [129-131]. 

Andreev sidebands in a system with two superconducting leads coupled by a processing 
spin proved important to study spin-transfer torques acting on the processing spin [132]. 
Spin-polarized Shapiro steps were studied in [133]. 


4. Andreev spectroscopy in F-S and F-S-F’ structures 

Andreev point contact spectra in S-F structures are modified with respect to those in S-N 
structures due to spin-filtering effects and the spin-sensitivity of Andreev scatfering [88,134-139]. 
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Spin-dependent phase shifts also crucially affect Andreev point contact spectra [72,140-145]. Point 
contacts have lateral dimensions much smaller than the superconducting coherence lengths of the 
materials on either side of the contact. Typically, a voltage will be applied over the contact, which 
makes it necessary to study in addition to the coherence amplitudes 7 and 7 also distribution 
functions. I follow the definition in Ref. [45,105], where 2x2 distribution function spin-matrices x 
and X for particles and holes are introduced which obey a transport equation 

\h{vp-V + dt)x-{-ioA + E)ox + xo{'yoA + S)^ (4.1) 

ih{vF- V-dt)x- {^oA + E)ox + xo{^oA + E)^ ^j-coW ^^ 2 ) 

The distribution functions matrices are hermitian, x = x^ and x = x^. The right hand sides of 
equation (4.1) and (4.2) contain collision terms (see Ref. [105] for details), which vanish in 
ballistic structures. In general, these distribution functions can be related to quasiclassical Keldysh 
propatators = —27ri(V oa;oV^ — J- 0 x 0 J^^) and = —27ri(V 0 x 0 — E o x oV^). The Fermi 
distribution function for particles, /p, and holes, /h, is related to x in the normal state by 
/p = (1 - x )/2 and /h = (1 - x)/ 2 . 

(a) Andreev processes in point contact geometry 

In this section I consider point contacts of dimensions large compared to the Fermi wavelength 
and small compared to the superconducting coherence lengths. In this case, the wavevector ky 
parallel to the contact interface is approximately conserved. The current on the ferromagnetic 
side of a point contact, being directed along the interface normal, can be decomposed into 

/ = /i-/r, + /ar (4.3) 


where the various terms are the incoming current, Jj, the normally reflected part, Jr, and the 
Andreev reflected part. Jar- The sign convention here is such that a positive current denotes 
a current into the superconductor. Thus, in the normal state the current I is positive when the 
voltage in the ferromagnet is positive. The various currents can be expressed as 


/x 




2Tih 


AcF 


d"g(k||) 

(27r)2 


de . 
yejx, 


(4.4) 


where e = — |e| is the charge of the electron, A is the contact area, and TIcf is the projection of the 
Fermi surfaces in the ferromagnet on the contact plane. For each value ky there will be a number 
of (spin-polarized) Fermi surface sheets involved in the interface scattering (in the simplest case 
spin-up and spin-down, or only spin-up), and the dimension and structure of the scattering 
matrix will depend on how many Fermi surface sheets are involved. The sum over a and 13 
runs over those Fermi surface sheets 1,..., v for each given value of ky. In the superconductor 
I assume for simplicity that only one Fermi surface sheet is involved for each ky. The reflection 
and transmission amplitudes for each ky are related to the scattering matrix as 

-Z) <“■=> 

where directions 1 and 2 refer to the superconductor, and 3 and 4 to the ferromagnet. R 12 is a 2x2 
spin matrix, R 34 is a 1 / x u matrix with elements Rap, T 14 is a 2 x r matrix with elements Tip, 
and T 32 is a R X 2 matrix with elements Ta2- The spectral current densities j'x are given by 

jR = \Ral3 - Ta2 «2 72 ^21 7l TipfSxp (4.6) 

/3 aP 

JAR = 1^02 V2 72 T2^^5xa, l>2 = (1 “ 72 .R 2 I 7l -R 12 ) 


OlO. 


(4.7) 
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Figure 6. (a) Shape function of the (spin-averaged) interface barrier potential. A shape parameter ct = 0... 0.7 Ap 
increases for increasing smoothness, (b) Spin-mixing angie i?s as a function of |k|| | for the shape functions in (a). i?s 
increases with increasing a. (c) The differential conductance of an F-FI-S structure with various degrees of smoothness 
if the FI barrier, at T = 0; cr increases from back to front in steps of 0.1 Ap. (d)-(e) Temperature dependence of the zero- 
voltage conductance of an FIM-FI-S point contact as predicted by (d) the modified BTK model [137] and (e) the spin-active 
interface model [82,140]. Adapted from [72,140]. Copyright (2010) by the American Physical Society. 


where 6 xp and 5xp are the differences between the distribution functions in the ferromagnet and 
in the superconductor. If there is no spin-accummulation present, they are independent on the 
index /3 and given by 

Sx{V, T; e) = tanh - tanh T; e) = Sx{V, T; -s), (4.8) 

with Tg the temperature in the superconductor. Equations (4.6)-(4.7) are valid for general normal- 
state scattering matrices S, and can be applied to non-collinear magnetic structures. For the 
case that all reflection and transmission amplitudes are spin-diagonal, and assuming on the 
superconducting side of the interface a spin-mixing angle i?, fhese expressions are explicitely 
given by 

jR = [I^^o-fI^ kt - + ko-P ki - Tte“‘’’7oP] Sx (4.9) 

ji = 2Sx, iAR = (k^4)^l7ol^ [ko-hl^ + ko-l^] <55 (4.10) 

with 7 o = —Aj{e -F if?), f? = -s/, wo± = (1 “ energy e is assumed 

to have an infinitesimally small positive imaginary part. Andreev resonances arise for energies 
fulfilling 1 = 7 o(£)^r^r^e^*’^ (in agreement with the discussion in section 2(b) for d — 0). 

On the other hand, for half-metallic ferromagnets the Andreev reflection contribution is zero 
in collinear magnetic structures. In non-collinear structures, however, the process of spin-flip 
Andreev reflection takes place, introduced in reference [140], and illusfrated there in figure 11. 
Spin-flip Andreev reflecfion is the only process providing particle-hole coherence in a half- 
metallic ferromagnet. Such structures are described by the theory deceloped in appendix C of [82]. 
Application of this theory to experiment on Cr02 is provided in [72,144]. A generalization to 
strongly spin-polarized ferromagnets with two itinerant bands is given in [140] with application 
to experiment in [141]. 

In figure 6 selected results are shown. In (a) and (b) it is demonstrated that the spin-mixing 
angle can acquire large values if a smooth spatial interface profile is used instead of an atomically 
clean interface. Correspondingly, in (c) Andreev resonances are more pronounced for smoother 
interfaces. In (d) and (e) a comparision of fhe model by Mazin et al. [137] for various spin 
polarizations P with the spin-mixing model for P=100% and various spin-mixing angles d shows 
that the two can be experimentally differentiated by studying the low-temperature behavior [72]. 

In an experiment by Visani et al. [146] geometric resonances (Tomasch resonances and Rowell- 
McMillan resonances) in the conductance across a Lao, 7 Cao. 3 Mn 30 /YBa 2 Cu 307 interface were 
studied, demonstrating long-range propagation of superconducfing correlations across the half 
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Figure 7. Illustration of notation used in text. For brevity of notation I sometimes omit the label 3 and 4, implying that a 
then means 3a and 0 means 4/3. E.g. in the right picture the ferromagnet has two spin Fermi surfaces (red and blue), 
labeled by a e {31, 3 4,} and /3 e {41,4},} etc. The superconductor's Fermi surface is drawn in green. 


metal Lao.yCao.aMnaO. The effect is interpreted in terms of spin-flip Andreev reflection (or, as 
named by the authors of [146], "equal-spin Andreev reflection"). 

Spin-dependent scattering phases qualitatively affect the zero- and finite-frequency current 
noice in S-F point contacts [73,74]. It was found that for weak transparency noise steps appear at 
frequencies or voltages determined directly by the spin dependence of scattering phase shifts. 


(b) Andreev bound states in non-local geometry 

A particular interesting case is that of two F-S point contacts separated by a distance L of the 
order of fhe superconducting coherence length. This is effectively an F-S-F' system, or if barriers 
are included, an F-Fl-S-FT-F' system. In this case, for a ballistic superconductor, one must consider 
separately trajectories connecting the two contacts [75]. Along these trajectories the distribution 
function is out of equilibrium, and equations (4.1)-(4.2) must be solved. In addition, the coherence 
functions at these trajectories experience both ferromagnetic contacts, and are consequenctly 
different from the homogeneous solutions 70 , 70 of all other quasiparticle trajectories. 

The current on the ferromagnetic side of one particular interface (positive in direction of the 
superconductor), can be decomposed in an exact way, 

I = h-^K +I'EO + (4.11) 

where the various terms are the incoming current, Ii, the normally reflected part, Jr, the Andreev 
reflected part, FaR/ arid the two non-local contributions due to elastic co-tunneling, /gc and 
crossed Andreev reflection, /car- 

There will be contributions from trajectories in the superconductor which do not connect the 
two contacts. These will be described by equations (4.6)-(4.7) above. Here, 1 will concentrate 
on the non-local contributions, which arise from the particular trajectories connecting the two 
contacts. Assuming the area of each contacf much smaller than the superconducting coherence 
lenght (however larger than the Fermi wavelength, such that the momentum component parallel 
to the contact interfaces are approximately conserved), one can identify all trajecfories connecting 
the two contacts, treating only one and scaling the result with the contact area. The solid angle 
from a point at the first contact to the area A' of the second contact is given by 517 = where 

A '2 is the projection of the area of fhe second contact to the plane perpendicular to the line 2, 2 ' 
which connects the contacts (see figure 7 for the notation). 

Using the conservation of ky, and that consequently d^S(k||) = d^S'(pF 2 )|n • vf 2 |/|vf 2 | = 
d^S'(pFQ,)|h • vfqI/IvfqI (where h is the contact surface normal), one can express fhe currents 
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Ix 


oo 

d^Si A2A'2 de . 

dfilpFa (27rfi)3L2 yJx, 


(4.12) 


where pf 2 is the particular Fermi momentum in the superconductor corresponding to a Fermi 
velocity in direction of the line 2,2' (I assume for simplicity that only one such Fermi momentum 
exists), and d^5'/df2 is the differential fraction of the Fermi surface of fhe superconductor per 
solid angle 17 in direction of the Fermi velocity vf 2 that connects the two contacts. Note that 
d^S/dil is the same at both contacts for superconductors with inversion symmetry, as then this 
quantity is equal at pf 2 and —pf 2 - Reversed directions are denoted by an overline: Pf 2 = —PF 2 
etc. Let us introduce scattering matrices S(12; 34) as well as S(21,43) at the left interface, the latter 
being equal to S(12; 34) for materials with centrosymmetric symmetry groups, which I consider 
here. Analogously, for the right interface let us introduce the scattering matrices S'(2'l'; 4'3') = 
S'(I'2'; 3'4'). The scattering matrices for holes are related to the scattering matrices for particles 
by S(21; 43) = S(21, 43)* etc. One obtains for this case 

+ Sxa (4.13) 

p a 

JR = X] l-^a/3 - Ta2 t'2 72 -R 2 I 7l Tip\^5xi} + ^ 7j Ri2 72 T2a\^5xa (4.14) 

a/3 pa 

JAR = V2 72 f2a\^5xa + ^ iTgi 7l Tip\^^Xp (4.15) 

““ ^ ^ 

JEC = V 2 U 22 ' T2,a'\^5Xa> + ^ iTgi Vj 7 j R 12 U 22 ' Rv Ti,p,\^5xp, (4.16) 

aa' pp' 

JCAR = X] V2 U 22 ' R2'1' 7i' f[,p,\^5xp> + ^ \Tpi vj 7j .Rj2 U 22 ' Ti'a'f^xa> (4.17) 

aP' pa' 


where the vertex corrections due to multiple Andreev processes are W 2 = (1 — 72 ^21 7i 1 ^ 12 ) ^ 
and nj = (1 — 72 ^2i)~^- ^or unitary order parameters {AA ~ (jq), let us define 17erg = 

[—— (e + iO^)^CTo]^ as well as 7 = —Al{e + il7), 7 = A/{e + il2). Using the amplitudes 


72' — 7?2'1' 7i' -Ri' 2': A' —^2'!' 7l' 77 i'2'i 
and denoting with L the distance between 2 and 2', 


U 22 ' = 


. — e 

-1 p 

C2' + 1 o ®2' 

J/2' 

, 72 = U22' 

. r^/ z^9/ “h E 

-1 r 

h'-' 172' 

, 72 = M 22 ' 


. A2' + /2'e 


172 


-S2' 


77/02' ~ i 


. ^ 2 ’ — R^ 2 ’^ 


17n, 


S2' 


(4.18) 

(4.19) 

(4.20) 


with C 2 ' = cosh(172'7//fit)F.2') and S 2 ' = sinh(172'L/fiiiF,2')- For f^e distribution functions one 
obtains for the two leads 


5xp = xa=3x{V,T;£), Sxa = 3xp = Sx{V,T;—e), (4.21) 

3 x^1 = 5xpi = x{V',T';£), 3xpi = 5xa'= 3x{V',T'-—e). (4.22) 

Here, Tg is the temperature in the superconductor, T and V are temperature and voltage in 
the left lead, and T' and V' are temperature and voltage in the right lead. The voltages are 
measured with respect to the superconductor. The expressions appearing in equations (4.13)- 
(4.17) have an intuitive interpretation, and selected processes are illustated in figure 8. These 
terms involve propagation of particles or holes, represented as full lines and dashed lines in 
the figure. Certain processes involve conversions between particles and holes, accompanied by 
the creation or destruction of a Cooper pair (loops in the figure), and correspond to the factors 
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Figure 8. Illustration of selected processes contributing to the expressions (4.13)-(4.17). Andreev reflections are 
denoted as loops, turning particles (p, full lines) into holes (h, dashed lines) or vice versa, (a)-(b) contributions to the 
reflection components, equation (4.14); (o)-(d) contributions to the Andreev reflection components, equation (4.15); (e)-(f) 
contributions to the coherent electron transfer components, equation (4.16); (g)-(h) contributions to the crossed Andreev 
reflection components, equation (4.17). 


71/ 7i'/ 71 / arid 7 j/ in (4.14)-(4.17). Propagation of particles or holes between the left and right 
interface is represented in these equations by the factors U 22 ' and U 22 ' ■ Vertex corrections V 2 and 
vi correspond to multiple Andreev reflections at either interface. The factors 72 and 72 combine 
propagation between the two interfaces with Andreev reflections at the other interface. 

As an example, for an isofropic singlet superconductor and collinear arrangement of the 
magnetization directions, one obtains 

jl = 4:Sx, jR,= |^2|i;+|^ |r^ - r^e"’707+P + 2|u-|^ |r^ - r^e“'’’707-l^] < 52 ; (4.23) 

jAR = mtif (|7+|^ + |7o|^) + (|7-l^ + l7ol^)] Sx (4.24) 

JEG = [mt\f\v+u+\‘^ {1 + l7ol‘‘(t'7i-l)^} + |l + |7ol'‘(i-t4)^}] (4-25) 

icAR = I 70 P [{ttt[f\v+u+\'^ {(4)^ + (^ 4 )^} + (*44)^1^-’^-!^ {(4)^ + (4.26) 

where 1 defined (with s = S2' and c = C2>) 

r± = 44e±‘’^'70, 7± = M± [r±c + is{A + r^e)/Q] (4.27) 

U±=[c-is{£ +r±A)in] 4 W± = [1 - 7±7o’"t’'4®"'"*’^r^' (4-28) 


Early studies of nonlocal transport in F-S-F structures include Ref. [76]. In Ref. [77] the non¬ 
local conductance was explained in terms of fhe processes discussed above for an F-S-F sfructure 
with strong spin-polarization. Andreev bound states appear on both ferromagnet-superconductor 
interfaces, which decay trough the superconductor towards the opposite contact. Parallel and 
antiparallel alignment of the magnetizations lead to qualitatively different Andreev spectra. The 
non-local processes have a natural explanation in terms of overlapping spin-polarized Andreev 
states. The density of states for the trajectory connecting the two contacts is obtained from 


N^{.e,x) = ^Re 


1 -7-t(£,a:)7-(-£* 
l-4 7+(£,a;)7* (-£* 


L — x) 
L — x)’ 


(4.29) 


and for the same expression holds with + and — interchanged. In figure 9 1 show an example 
of such a sefup. As can be seen, avoided crossings of Andreev boimd sfates with equal spin 
polarization play an important role in such systems. At the avoided crossing the bound states 
have equal weight at both interfaces. This contrasts the case when bound states have opposite 
spin polarization, where no avoided crossings appear, and the case when the two S-F interfaces 
have markedly different spin-dependent phase shifts, in which case bound states do not overlap 
and stay localized at one of fhe two interfaces only. The spin-polarization and weight of the 
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Figure 9. Andreev bound states in an F-S-F structure of the type shown in figure 7. (a): P = as function of e and 

at a position in S midway between the contacts, for L = 2go (with the coherence iength of the superconductor = 
and i?L = 0.77r. An avoided crossing appears for equaliy spin-polarized Andreev states, which is absent for 
opposite polarization, (b) dependence on L/^o for fixed i?r = = O.Ttt. (c)-(f): P as function of e and xior L = 2^0 

and = O.Ttt, and (c) »9 r = O.Ttt (d) i7r = O.Stt (e) i?r = —O.bTr (f) i9r = —O.Ztt. At the avoided crossing all bound 
states have equal weight at both interfaces. In all other oases the bound states for fixed spin projection are localized at 
one interface only. In all panels = (r’^r|)R = 0.95. 


bound states at energies Sb arid —Sb determine the magitude of the nonlocal currents due to 
crossed Andreev reflection and elastic co-tunneling. A detailed discussion of how the weights, 
transmission probabilities, and bound-state geometries influence CAR and EC processes is given 
in [77]. In diffusive S-F systems a theory for nonlocal transport was developed in [147]. 

Equations (4.23)-(4.26) have been applied to the study of thermoelectric effects in non-local 
setups [148,149]. The contributions to the energy current are obtained as 

le = 4i - 4 r + 4ar - 4ec + 4car, (4.30) 

where the respective contributions are given by analogous equations as in equations (4.12)-(4.17), 
however with replacing the charge e by energy e. Spin-dependent interface scattering phases in 
combination with spin filtering leads to giant thermoelectric effects in F-S-F devices [148,150]. 

There have been a number of experimental studies of non-local transport in S-F hybrid 
structures, e.g. [58,151,152]. Cold et al. investigate S-FF-S junctions with two parallel running 
wires [153]. A very recent experiment shows a non-local inverse proximity effect in an (F-N-F')-S- 
N structure [154]. The inverse proximity effect transfers magnetization from the ferromagnet into 
the superconductor or across a superconductor. For strongly spin-polarized systems this occurs as 
a result of spin-mixing phases [89,107]. A combination of non-local effects in S-F structures with 
non-equilibrium Andreev interferometer geometries, in analogy to experiments in S-N structures 
[155,156], seems to be another exciting avenue for future applications [157]. 


5. Generalized Andreev Equations 

In this section 1 discuss the physical interpretation of the coherence functions. To this end, 1 
present a generalized set of Andreev equations which is equivalent to equations (3.9)-(3.10). Let 
us define for each pair of Fermi momenta —pp and corresponding Fermi velocities vp{pp), 
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R{p) = Ro+hVF{PF){p-Po), Rip) = Rl-fiVF{PF){p-pi), P 0 <P<P 1 - (5-1) 

Using dp = hvF- V, let us define the following differential operators 


V = 


-\dp + r 
-A 


A 

idp-S 


V = 


-idp + r 

-z\ 


idp-E 


(5.2) 


which fulfill T) = —txDti. Let us also define fhe adjoint operator 'D*{p, dp) = {p, —dp). For a 

fixed conjugated trajectory pair the set of generalized Andreev equations (retarded and advanced) 


IS, 




V 


w'^(Pi) = -71 °w”(pi) 
w"(Po) = -7 o°m”(Po) 

v'^{po) = -^l°u^{po) 


(5.3) 


(5.4) 


where the boundary conditions at p = po and p — pi for the solutions fulfill the restrictions shown 
on the right hand side of the equations. Then the relation between the Andreev amplitudes u, v, 
u, and V and the coherence amplitudes 7 and 7 is given along the entire trajectories by [105] 


~R,A R.A ■’R^A R,A ~ R,A R,A 

V = —7 ou , V = —7 ou 


(5.5) 


with 7 ^ = 7 , 7 ^ = 7 , 7 '^ = 7 ^, 7 '^ = 7 ^. It is easy to show that the following conservation law along 
the trajectory holds 


rs o 


= 0 . 


(5.6) 


Thus, the matrix inside the curly brackets is given by its value at one point on the trajectory. 
The off-diagonal elements are zero due to the conditions in Eq. (5.3) for po and pi, leading 
fo (M^)f ov^ = (ti^)f ou^ and ov^ = (0'^)^ ou” along the entire trajectory. The diagonal 

components dp[(u'^)f o u” — (u'^)f o u”] = 0 and dp[(u^)f o u” — (u'^jf o u”] = 0 translate into 
dp[(u'^)f o (1 — 7” o 7”) o u”] = 0 and dp[{u^)^ 0(1 — 7” ° 7 ”) ° — 0. In particular, for Andreev 

bound sfates 1 — 7” o 7” = 0 or 1 — 7” o 7” = 0 , and thus this property is conserved along the 
entire trajectory. 

If one writes Eq. (5.3) formally as 27 o (7 = eU, then the conjugated equation V oU = —ell holds 
with U = tiUti, which leads, however, to a system identical to Eq. (5.3). The adjoint equation 
27* o U = eU defines adjoined Andreev amplitudes (left eigenvectors) u, v, u, and u. These are, 
however, equivalent to the advanced eigenvectors in Eq. (5.4). 


6. Conclusion 

I have presented theoretical tools for studying Andreev reflection phenomena and Andreev 
bound states in superconductor-ferromagnet hybrid structures. Concentrating on ballistic 
heterostructures with strong spin-polarization, I have formulated theories for point contact 
spectroscopy and for nonlocal transport, as well as for Andreev states in Josephson structures in 
terms of coherence functions and distribution functions. The connection to coherence amplitudes 
appearing in the solutions of Andreev equations has been made explicit. The formulas for non¬ 
local transport have been given in a general form, allowing for non-collinear geometries, and 
using the normal-state scattering matrix as input. 

Data accessibility. All data can be obtained from the author. 
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